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ABSTRACT 

Pitch-angle scattering rates for cosmic-ray particles in magnetohydrodynamic (MHD) simulations 
with imbalanced turbulence are calculated for fully evolving electromagnetic turbulence. We compare 
with theoretical predictions derived from the quasilinear theory of cosmic-ray diffusion for an ideal¬ 
ized slab spectrum and demonstrate how cross helicity affects the shape of the pitch-angle diffusion 
coefficient. Additional simulations in evolving magnetic fields or static field configurations provide 
evidence that the scattering anisotropy in imbalanced turbulence is not primarily due to coherence 
with propagating Alfven waves, but an effect of the spatial structure of electric fields in cross-helical 
MHD turbulence. 

Subject headings: diffusion — scattering — magnetohydrodynamics (MHD) — (ISM:) cosmic rays 


1. INTRODUCTION 

Both turbulence and magnetic fields are phenomena 
which permeate essentially every subdiscipline of astro¬ 
physics. Studying the combined effect of both phenom¬ 
ena on fundamental processes such as the diffusion of 
charged particles in the context of magnetohydrodynam¬ 
ics is therefore a task of utmost importance for gaining 
a better understanding of the universe. 

The magnetohydrodynamic (MHD) turbulence of 
many astrophysical systems exhibits high degrees of cross 
helicity, which means that fluctuations of the bulk ve¬ 
locity of the medium and of magnetic fields induced in 
the plasma are strongly correlated. This property, which 
is also known as Alfvenicity w hen referring to the solar 
wind (Matthaeus et al. 2004) or as i mbalanced turbu¬ 
lence in the astrophysical coirimunity (Perez & Boldyrev 
2009 1, has been measured in observations of the fast solar 
wind (Marsch & Tu 1990) and is i mportant upstream of 
super nova remnant (SJNR ) shocks (Schlickeiser & Shalchi 
2008) as well as a driving force of t he dynamo effect in 
accretion disks and young galaxie s (Yoshizawa & Yokoi 
1993 Brandenburg fc Urpin||1998 ). 


The rate at which MHD turbulence scatters the pitch- 
angle of charged particles is a crucial parameter in many 
of these settings, particularly for diffusive shock acceler¬ 
ation in SNR shock fronts or measurements of velocity 
statistics in the solar wind. Eor example, the isotropiza- 
tion of cosmic-ray particles at SNR shock fronts is often 
considered to occur at the rate of Bohm diffusion, which 
can o nly be strictly derived in the limit of strong turbu¬ 
lence (Shalchi 2009). In weak or intermediate magnetic 
turbulence with dn < Bq, isotropization will be signif¬ 


icantly slower and the efficiency of the diffusive shock 
acceleration process will be diminished. 

Pitch-angle scattering on short timescales is particu¬ 
larly impo rtant in quasipe rpendicular shocks. It has been 
proposed (jSagdeev 1966) that low-energy ions entering 
the shock toot troni the upstream region can be trapped 
by the electrostatic potential drop in the ramp. Within 
a few gyroperiods, ions may be accelerated while surfing 
along the convective electric field parallel to the surface 
of the shock front. Thus they can gain even higher en¬ 
ergies, which are required for subsequent efficient Eermi 
acceleration. The details of this process, however, de¬ 
pend on how fast pit ch-angle scattering pro ceeds over 
these few gyroperiods ( Kirk fc Heavens|[l989 ). 

The acceleration of dust particles moving through the 
interstellar medium at velocities close to the Alfven speed 
similarly depends on the rate of pitch- angle scattering 
(Lazarian & Yan|2002 Hoang et al. 2012). Consequently, 
a large number or comparisons between analytical predic¬ 
tions of thi s rate and test-particle simulations has been 
published ( Michalek & Ostrowski 1998 ^in fc Shalchi] 


2009 fDalena et ai.||2012' 'i'autz et al. 112013” 

Most of these studies assumed a power-law spectrum 
for the magnetic turbulence and neglected the influence 
of cross helicity. Starting with a random distribution of 
the complex phases of the individual wave modes, real¬ 
izations of the magnetic field were generated from this 
analytic description and used to compute test-particle 
trajectories. While this scheme represents a computa¬ 
tionally efficient method of creating turbulence with a 
wide inertial range, it is only useful for comparisons with 
a magnetostatic turbulence model as it fails to capture 
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the nonlinear effects important for a realistic description 
of the dynamic evolution of the turbulence, an d in par¬ 
ticular of the electric fields (Arzner et al 


1 [ 200 ^ . 

using fie 


These shortcomings can be avoided by using helds gen¬ 
erated in numerical three-dimensional MHD simulations 


to compute the test-par ticle scattering (e. g. Lehe et al. 


2009 Lynn et al. 2012|), at the cost of a narrower iner- 
tiai range than is possible with the analytic approach. If 
one aims to investigate the effects of cross helicity, which 
manifest themselves primarily in the dynamic propaga¬ 
tion of Alfven waves and a reduction of the motional elec¬ 
tric field Binot = — u X b, a full-MHD approach becomes 
necessary to obtain realistic results at particle speeds 
comparable to the Alfven velocity. 

Only few authors have numerically investigated the 
scattering of cosmic-ray particles in plasmas o f non¬ 
vanishing cross helicity. Beresnyak et al. (2011) stud¬ 
ied the scattering of relativistic test-particles m balanced 
and imbalanced MHD fields, but did not report observ¬ 
ing any influence of the cross h elicity on the pitch -angle- 
averaged scattering frequency. Teaca et al. (2014) found 
that the decreased magnitude of emot in cross-helical 
MHD turbulence slows down the isotropization of Alfven- 
speed particles significantly. 

Here we present our results on the pitch-angle depen¬ 
dence of test-particle scattering in simulations of cross¬ 
helical MHD turbulence and compare them to predic¬ 
tions of the pitch-a ngle diffusion coef ficient derived from 
quasilinear theory ( Schlickeiser||1989 ). 

This article is structured as follows: After an overview 
of various turbulence models used in quasilinear theory 
and how they affect the predicted pitch-angle diffusion 
coefficient in imbalanced turbulence (Section]^, we dis¬ 
cuss the computational methods employed in our inves¬ 
tigation (Section]^, including a detailed analysis of the 
structure of electric and magnetic fields in MHD simula¬ 
tions of plasma turbulence with various cross helicities. 
Having explained the rationale behind our choice of diag¬ 
nostics, we present the results on pitch-angle scattering 
for test-particle simulations in these cross-helical turbu¬ 
lent fields (Section]^ before we summarize our findings in 
order to draw several important conclusions (Section]^. 


2. QUASILINEAR PREDICTIONS 
2.1. Magnetostatic turbulence 


The diffusion of a beam of charged cosmic-ray parti¬ 
cles with velocity u in a turbulent magnetic held b(x) is 
commonly described in terms of a Fokker-Planck equa¬ 
tion for the cosmic-ray phase-space density /(/i, z), where 
fi = cos(u^/'!;) is the cosine of the pitch-angle and z is 
the c oordinate along the global mean magnetic held Bqz 
(e. g. Schlickeisei]|2002 Slialchi|[2009 1: 


dt 


+ Vfl 


dz 


A 

dfi 



( 1 ) 


The form of the quasilinear pitch-angle diffusion coef- 
hcient = limt_>,oo — ir(0)]^)/2t was hrst derived 
for cosmic- r ay par ticles in magnetostatic slab turbulence 
by Jokipii (|1966 ). Expanding the turbulent magnetic 
held b(xj = Bqz + db(z) in a perpendicular perturba¬ 
tion 5bj_(z) T z that only depends on the coordinate 


along the zeroth-order mean held Hgz, he found that 

D,, = ^ D2-(l - ^2) ^ - l)fc-i. (2) 

Here, the power spectrum of the magnetic held pertur¬ 
bations is assumed as a slab spectrum \5hj[^{kz)\^Akz = 
^B^(s — ^ spectral index s, a lower 

cutoh wavenumber at fcmin, and root-mean-square am¬ 
plitude ((5b^)^/^ = |5B|. 

As the typical velocity of cosmic-ray particles u ~ c 
exceeds the typical Alfven velocity in the interstellar 
medium va ^ lOkms”^ by several orders of magnitude, 
the magnetic turbulence can be approximated as being 
time-independent. In the limit of very small wave fre¬ 
quencies D ^ u/cniin, we Can also presume that the char¬ 
acteristic electric held strength 6E becomes negligibly 
small since Faraday’s law implies SE ^ D|dB|/fc. As a 
consequence of these assumptions, the cross helicity does 
not enter into equation (§. 


2.2. Electrodynamic turbulence 

For charged particles with v ^ va-, however, the prop¬ 
agation speed of shear-A lfven waves must be consid¬ 
ered. As hrst shown by Schlickeiser (1989), all pow¬ 
ers of 7 = va/v must be included in the derivation 
of to obtain the correct result. Accounting for 
the time-dependence of the magnetic turbulence and, 
through Faraday’s law, the resulting electric helds self- 
consistently, one hnds for the quasilinear pitch-angle dif¬ 
fusion coefficient in incompressible MHD turbulence: 


Dfifi — ^ ^ Bq (1 /i ) 

/3=±1 


Here /? S {-1-1,—1} distinguishes shear-Alfven waves 
traveling along (/3 = —1) or opposite (/? = -1-1) the mean- 
held axis, the only types of waves that we consider in 
our simplihed model of incompressible slab turbulence. 
We neglect additional contributions to the Fokker-Planck 
equation proportional to dpf and d^f that result in adi¬ 


abatic focusing and stochastic heating, focusing only on 
the pitch-angle evolution of a monoenergetic particle dis¬ 
tribution. 

The total turbulence energy of the system can be sep¬ 
arated into the energies contained in counter- and co¬ 
propagating waves traveling opposite to or along the di¬ 
rection of the mean magnetic held, the so-called positive 
and negative Elsasser energies = ((u ± (5b)^)/4. In 
quasilinear theory, an analogous decomposition into two 
Alfven-wave populations is used for only the magneti c 
turbulence energy as well (|Dung & Schlickeiser||1990|), 
and we have written the respective components of the 


magnetic-turbulence power spectrum as £g{kz) above. 
The wave numbers resonant with each wave population 


are dehned as k^es = + /3va)- In case the power 

spectra for both wave populations are not equal, the de¬ 
gree of imbalance between both propagation directions 
can be captured in the normalized cross helicity 




£+ -£- 
£++£-' 


( 4 ) 
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Written in terms of CTc, in the electrodynamic case the 
pitch-angle diffus ion coefficient becomes (Dung & 
Schlickeisei]|1990 ) 

ji ^ ^ ^ n w I I • I ^ 

D„ 




3=± 


s'3-l 

min 


X v" {l + P^J.'J) |/r-k/37|® {1+Pac). (5) 


To obtain this equation, we assume that the magnetic 
turbulence spectrum can be decomposed into contribu¬ 
tions by CO- and counter-propagating Alfven waves in 
the same manner as the Elsasser energies, so that the 
two spectra 

|(ib±(fc,)|2dfc, = {SB^f (1 - S±) Ktin" dfc, (6) 

describe the magnetic turbulence due to the two Alfven- 
wave populations, where = 5B^(1 ± tJc) is the 

root-mean-square of the magnetic turbulence in the two 
wave populations. 


2.3. Magnetodynamic turbulence 

In order to distinguish between the pitch-angle scat¬ 
tering caused by the electric-field acceleration and that 
due to the magnetic field alone, we have performed test- 
particle runs in which the electric-field component of the 
Lorentz force was not included. Without electric heat¬ 
ing, the pitch-angle diffusion equation (f^ becomes exact 
again, although must be adapteoto this ‘magne- 
tody namic’ turbulence m odel. Generalizing the results 
that Shalchi et al. (2009) published for balanced turbu¬ 
lence td''arBirtrary cross helicity, we obtain for the pitch- 
angle diffusion coefficient: 


/3=±1 


n 








X {1 — g?) l/i-h/Syl® ^ {l + jdac)- (7) 

As one can easily check, eqns. ([^ and ^ are identi¬ 
cal to (H) if the Alfven velocity becomes n^igibly small 
compared to the particle velocity {v ^ va or equivalently 
7 ^ 1 ) and both wave populations exhibit the same spec¬ 
tral exponent (s’*" = s~). The quasilinear predictions for 
the pitch-angle diffusion coefficient in all three different 
cases in zero cross-helicity turbulence and strongly im¬ 
balanced turbulence are compared in Figure [l] 


3. METHODS 
3.1. MHD setup 


In order to compare the validity of the quasilinear slab- 
turbulence results for the pitch-angle diffusion with nu¬ 
merical r esults, we use the pseudospectral MHD code 
Turbo (Teaca et al. 2009) to obtain representations 
of turbulent fields in which we propagate test-particles. 
Employing Alfvenic units. Turbo solves the equations 
of resistive incompressible MHD on a cubic grid with pe¬ 
riodic boundaries: 


dtu = -(u • V)u -I- (b • V) b -K u V^u-fP - Vp, (8) 
a4b = -(u-V)b + (b-V)u-Kr7V2b+f'' (9) 


Here p = p/p is a re-normalized pressure that ensures 
that the velocity field remains divergence-free, v and 




Fig. 1. — Shapes of the pitch-angle diffusion coefficient ac¬ 
cording to quasilinear theory for Alfven-speed particles {v = va) 
in magnetostatic turbulence (dotted), magnetodynamic turbulence 
(dashed), and fully electrodynamic MHD turbulence (solid lines) 
with amplitude SBJBq = 0.1 and Kolmogorov-like spectra s+ = 
s~ =5/3 for normalized cross helicities CTc = 0.0 and Cc = 0.9. 
(A|| = 27r/fcniin) 

r] are kinetic viscosity and magnetic diffusivity, respec¬ 
tively, while f„ and ff, are forcing fields. The magnetic 
field b includes a constant guiding field Bq applied along 
the 0 direction. 

The cross-helical forcing scheme we use decomposes 
the Fourier-transformed fields u(k) and b(k) into eigen- 
modes of the curl operator. These helical eigenmodes 
for wave vectors in the ellipsoidal shell 2.5 < k'j_ + 
^ 3.5 are then forced separately in such a way 
that the injected amounts of energy (£), kinetic helic¬ 
ity ("Hkin), magnetic helicity {'Hma.g), and cross-helicity 
{JC = (u • Jb) = ac£) are constant. Accordingly, we 
define the perpendicular and parallel correlation lengths 
as A_l = 2n/3 and Ay = 10 A_l. The A:-space ellipsoid is 
shortened by a factor of 10 along the kz axis to allow for 
anisotropic forcing of the turbulent fields so that the tur¬ 
bulence can attain stable steady-states in the presence of 
a magnetic guide field. In particular, turbulent energy 
can be supplied to the system at the rate e = d£ jdt with¬ 
out affecting any of the helicities, or cross-helicity can be 
injected at a desired rate cr = ((?A/i9t)/e to drive the 
system towards an imbalanced steady-state with /C yf: 0. 

The fields used for the test-particle runs were obtained 
on 512^ grids with physical side lengths of 27r x 27r x 207r 
and periodic boundary conditions. This anisotropically 
shaped simulation domain allows the turbulence to reach 
a steady-state in the presence of a strong magnetic mean- 
field along the z direction, Bq = 10|(5B|, with |(5B| = 

((b-Bo)2)i/2. 

Viscosity and magnetic diffusivity are chosen to ensure 
that Kolmogorov’s dissipation length Ik = is 

resolved by the numerical grid in all directions. Since we 
use isotropic viscosity and diffusivity coefficients, energy 
dissipation remains an essentially isotropic process, while 
the physical size of the numerical domain is anisotropic. 
Kolmogorov’s dissipation length therefore corresponds 
to fewer grid cells in the direction along the mean-field 
than perpendicular to it, and resolving it with the same 
accuracy in both directions is impossible. With our 
choice of dissipation coefficients, Kolmogorov’s dissipa¬ 
tion length Ik fulfills 12 and ~ 1-2, 

where k^^ = tt/Ax and fcAy = tt/Az are the wavenum¬ 
bers corresponding to the grid resolution in perpendicu¬ 
lar and parallel direction, respectively. 

Driving turbulence with an energy-injection rate e = 
0.1 and cross-helicity injection rate cr = 0.0 leads to a bal- 
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Fig. 2.— Steady-state spectra of positive (black) and negative 
(red) Elsasser energies, = (u di (5b)^/4, for balanced (ctc = 0.0, 
dashed lines) and strongly imbalanced (tJc = 0.9, solid lines) cases, 
a) Spectrum perpendicular to Bq. b) Parallel spectrum 


anced steady-state with ac = JC/S — 0.0, while changing 
the cross-helicity injection to cr = 0.8 eventually yields a 
steady-state cross helicity of CTc = 0.9. The latter value is 
obtained from averaging the scalar product /C = (u • ^b) , 
which saturates at /C = 0.9 £ when the dissipation of 
energy and cross helicity matches their injection rates. 
Similarly, injecting positive or negative cross helicity at 
a lower rate {a = ±0.5) results in a steady-state config¬ 
uration with (Jc = ±0.6, respectively. 

Figure [^provides comparisons of steady-state spectra 
for Ctc = 0.0 and (Tc = 0.9. For the strongly imbalanced 
case (solid lines), the Elsasser energy with the same (in 
this case, positive) sign as the total cross-helicity dom¬ 
inates over the opposite Elsasser energy at low k < 20, 
while both Elsasser energies are in equipartition in the 
balanced case. The positive Elsasser energy is gener¬ 
ally associated with waves propagating opposite to a 
strong magnetic mean-field. We deduce that the major¬ 
ity of the propagating shear Alfven waves in the imbal¬ 
anced case has a phase velocity directed opposite to the 
mean field (their energy is respresented by while the 
wave distribution in the balanced case exhibits no such 
anisotropy. 

The spectra in the direction parallel to the mean field 
show a similar dominance of the Elsasser energy with 
the same sign as the cross-helicity (Figure]^. As ana¬ 
lytical predictio ns and recent s olar-wind measurements 
have shown (see Podesta (20091 and references therein), 
one expects a scaling £ oc k~^, and indeed this scal¬ 
ing fits our simulations well in the balanced case. In the 
cross-helical cases, the spectral indices of t he positive and 


negat ive Elsasser energies differ slightly (Grappin et al. 


19831. 

The imbalance of the cross-helical cases is also obvious 
in plots of the spatial distribution of the cosine between 
the velocity field and the magnetic field b = Bq ± ^b. 
The zero-cross-helicity case exhibits positive and nega¬ 
tive values of cosZ(u, b) in approximately equal propor¬ 


tion (Figure [^,b). For the strongly cross-helical case 
with (Tc = 0.97 the velocity field shows a similarly weak 
correlation with the total magnetic field (FigurejS:; how¬ 
ever, note the small changes in the slope of me his¬ 
togram curve at the extremal values cosZ(u,b) = ±1 in 
b), while positive values of cosZ(u, 5b) dominate if the 
alignment between the velocity field and only the fluctu¬ 
ations of the magnetic field is considered (Figure [^). A 
similar plot of cosZ(u, 5b) for the balanced case is not 
depicted since it varies only negligibly from the total-field 
alignment already shown in Figure!^. 

Due to the almost identical distribution of the align¬ 
ment of u and (Bq ± 5b) in balanced and cross-helical 
turbulence, both cases exhibit a similar level of isotropy 
of the motional electric field along the mean-field di¬ 
rection. The distribution of is symmetric around 
6^ = 0 with an approximately expo nential drop-off in 
both directions (Breech et al. 2003), with the global 
average (e^) = O.OOU458(|e|) in the balanced case and 
(cz) = —0.000034(|e|) in the strongly cross-helical case 
((7e = 0.9). 

In order to investigate the presence of incompressible 
Alfven waves, we switched off forcing and dissipation 
once the balanced and imbalanced field configurations 
had reached a steady-state. Setting f„ = ff, = ry = = 0, 

we ensure that the propagation of shear-Alfven waves 
is not disturbed by non-ideal effects, which are not ac¬ 
counted for in the derivation of the quasilinear diffu¬ 
sion coefficients. Comparing the power spectra of the 
magnetic field Fourier-transformed in space and time 
(Figure |4), we find evidence of Alfven waves traveling 
along and opposite the mean-field direction and there¬ 
fore confirm that the evolution of the MHD turbulence is 
shaped mainly by propagating shear-Alfven waves. The 
x-shaped dispersion plot of the balanced case is symmet¬ 
ric with respect to the k^ = 0 line, indicating that the 
waves propagating parallel and antiparallel to the mean 
field are equally strong. In the strongly cross-helical case 
(iJc = 0.9), there is more energy in the modes along 
u! = —VAkz than in the modes with the oppositely di¬ 
rected phase velocity along u) = +VAkz, particularly at 
larger kz, indicating once more a dominance of Alfven 
waves propagating opposite to the mean- field. 


3.2. Test-particle setup 

We continue using zero forcing and dissipation while 
we compute the forces that test-particles experience in 
these turbulent steady-state fields. The trajectories of 
these test-particles are evolved using an implicit fourth- 
order Runge-Kutta solver with an adaptively determined 
timestep size. Their velocity at each timestep is calcu¬ 
lated using the Lorentz force computed from the MHD 
magnetic field b and, in some runs, the motional electric 
field ejiiot = — u x b, as interpolated from the grid data 
using third-order spline functions. 

At the starting time t = 0 of the test-particle simula¬ 
tion, ensembles of Ap = 2 x 10'^ particles are distributed 
randomly on the grid, with an energy equal to Ep = 
v{0Y/2 and a velocity Vz = along the direction of 

the mean-field magnetic field. The initial pitch-angle co¬ 
sine fi differs for each ensemble, taking values from -1.0 
to 1.0. To demonstrate the interaction with propagating 
Alfven waves, we initialize the test-particle ensembles 
with velocities close to the Alfven velocity, using both 
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Fig. 3.— a) Alignment of velocity and magnetic fields, cosZ{u, b), in balanced turbulence, b) Histograms of alignment in steady-state 
balanced turbulence (a), alignment of velocity and total magnetic field (c) and of velocity and magnetic fluctuations (d) in strongly cross¬ 
helical turbulence, c) Alignment distribution of cosZ(u, b) in strongly cross-helical turbulence (gc = 0.9). d) Distribution of cosZ(u, (5b) 
in the same snapshot of strongly cross-helical turbulence 
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Fig. 4.— Power spectra P(kz,io) of a) balanced {gc = 0.0) and 
b) strongly cross-helical (gc = 0.9) MHD turbulence 

i;(0) = VA and n(0) = 3 va - The charge-to-mass ratio of 
the test-particles is chosen such that their initial gyrora- 
dius is rg{0) = v{0)/{qBo/m) = 0.12A_l for Alfven-speed 
particles and rg(0) = 0.3(3A_l for w(0) = 3va- 
Since our goal is to compare analytical pitch-angle dif¬ 
fusion coefficients predicted by quasilinear theory with 
simulation results, we define /r as the cosine of the an¬ 
gle subtended by the particle velocity and the direction 
of the global mean magnetic field (the 2 ; direction) and 
thus follow the QLT convention. Using an only locally 
averaged mean-field direction as reference would yield a 
pitch-angle diffusion coefficient that would not fully con¬ 


form to the quasilinear definition. QLT is derived by 
expanding the equations of motion in a small perturba¬ 
tion of a globally constant magnetic field, hence it is the 
global mean-field direction that is relevant for a compar¬ 
ison. 


3.3. Measuring pitch-angle scattering 

In order to compare the pitch-angle scattering of test- 
particles on short timescales with quasilinear estimates, 
we define 




Tg / 


( 10 ) 


as the scattering rate after one gyroperiod. Here /r(0) 
and /r(Tg) are the pitch-angle cosine at the test-particle 
injection time and after one gyroperiod, respectively. 
Whereas as defined above in the context of QLT, 
corresponds to an infinite-time limit, we decided to mea¬ 
sure instead for several reasons. Although calculat¬ 
ing Agg for timescales significantly longer than Tg could 
be expected to improve the agreement with quasilinear 
calculations of Dg_g, we will demonstrate that this is not 
the case, due to various effects which were not included 
in the derivation of the quasilinear coefficients above, but 
which are inevitably present both in our simulations and 
in realistic MHD turbulence. In Appendix]^ we demon¬ 
strate that Agg is bounded below by 
First, as the pitch-angle cosine fj, is limited to the inter¬ 
val [—1, -|-1], its evolution is subject to boundary effects. 
Measuring the mean-square pitch-angle deviation 

A^2(t) = ([Mt)-M0)f) (11) 

for large t will, in general, underestimate the diffusion 
coefficient Dgg{fi{0)) at the initial /r(0) as some test- 
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particles will have been reflected at the boundaries. Even 
before this reflection occurs, however, varies 

strongly with /r and is therefore difficult to determine 
if the pitch-angle distribution spreads too quickly. This 
is particularly relevant if the turbulence amplitude is 
too large (for instance 6B/Bq ^1). In that case, even 
particles with a gyroradius as small as rg{0) = 0.36A_l 
isotropize almost completely within a single gyroperiod. 
Therefore, we only present results for an intermediate 
turbulence strength of SBjBo — 0.1. 

Second, our goal is to investigate the role of electric 
fields, which accelerate or decelerate test-particles on 
longer timescales. Since the Fokker-Planck equation Q 
does not account for these effects, good agreement can 
only be expected if the timescales on which is mea¬ 
sured are short enough that the change in kinetic energy 
is negligible. 

Most importantly, we will demonstrate that pitch- 
angle scattering on longer timescales is severely con¬ 
strained by the approximate conservation of the magnetic 
moment, 

^2 


the fir st adiabatic invariant in slowly varying magnetic 
fields (Landau fc Lifshi'tz]|1982 ). The standard formula- 


proximate conservation of (Goldstein 

o 

00 , 

Jaekel 

1998 see however 

Schlickeiser & Jenko 

21 

JTDI. 

As an 


the evolution of on timescales exceeding several gy- 
roperiods, so that A^^, which is computed after only 
one gyroperiod, is not affected by these constraints and 
therefore allows for a better comparison with quasilin- 
ear predictions. Although the absolute value of A^^ will 
not necessarily be close to the value of in all cases, 
the scaling with the pitch-angle cosine fj,, which is our 
primary concern in these investigations of cross-helical 
turbulence, will be almost identical. 


4. PITCH-ANGLE SCATTERING IN CROSS-HELICAL 
TURBULENCE 

4.1. Fully electrodynamic runs 

We begin with the most realistic scenario, fully electro¬ 
dynamic turbulence, in which the evolution of the MHD 
fields is computed in parallel with the test-particle tra¬ 
jectories. The initial conditions corres pon d to the field 
configurations presented in Subsection |3.1[ In Figure 
we show the evolution of the kinetic energy, the mag¬ 
netic moment, and the mean-square pitch-angle devia¬ 
tion for two representative ensembles with v{0) = ivA 
and p,(0) = ±0.6 in balanced and strongly imbalanced 
turbulence. 

For both ensembles in the balanced-turbulence case, 
energy and the averaged magnetic moment oscillate with 
the gyrofrequency, but the values of the minima attained 
after each full gyroperiod show only negligible variation. 
The mean-square deviation of Mmag, as indicated by the 
shaded regions, converges to about 12 % of the initial 
value as the magnetic moment is conserved to a large 
degree because of the relatively weak perturbation am¬ 
plitude SB/Bq. This puts a limit to the evolution of 
however, since the kinetic energy, the magnetic 
moment, and the pitch-angle cosine are approximately 
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Fig. 5. — Evolution in fully electrodynamic turbulence for the 
relative change of the kinetic energy AE, the ensemble-averaged 
magnetic moment ALmag (Icr-deviation shaded in blue/red), and 
the mean-square displacement of the pitch-angle cosine, A/r^, for 
particles with v{0) = 3 va and /r(0) = —0.6 (red) and /r(0) = +0.6 
(blue) 
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(right), for four values of a 




















































7 



Fig. 7.— Fully electrodynamic case: pitch-angle scattering rate for cosmic-ray particles with v = and gyroradius ro = 

v/{aBo) = 0.12Ay determined numerically (solid dots) and best-fit curves from eqn. (5) (dotted lines), and for particles with v = Sva and 
gyroradius rp = 0.36Ay (hollow dots and dashed line), a) Balanced turbulence with = 0.0 b) Strongly cross-helical turbulence with 
(Tc = 0.9 c) Cross-helical turbulence with CTc = —0.6 d) Cross-helical turbulence with CTc = 0.6 


connected via 

E « + Efi^. (13) 

This relation is valid if we assume that the magnetic tur¬ 
bulence is so weak that the difference between the aver¬ 
aged velocity components (ujj) along the local magnetic- 

field direction and (u^) along the global mean-field can be 
neglected, the difference being of order {SB/Bq)'^. Con¬ 
sequently, the minima of after each gyroperiod 

increase in value for two gyroperiods before settling on 
an almost constant value. 

In imbalanced turbulence with (t^ = 0.9, the test- 
particle ensemble with a positive pitch-angle of ^(0) = 
-1-0.6 evolves similarly to the balanced case. However, 
the counter-propagating particles with /r(0) = —0.6 are 
much closer to resonance with the dominating Alfen-wave 
population, which is traveling opposite to the mean-field 
direction as well. Consequently, these particles undergo 
deceleration as a result of inverse Landau damping and 
their kinetic energy exhibits a steady decrease (red lines 
in Figure [^. While the magnetic moment is approxi¬ 
mately conserved again, the variation in energy makes it 
possible for pitch-angle scattering to continue on longer 
timescales; hence the minima of A/i^(t) after each gy¬ 
roperiod grow steadily. Obviously, a comparison of the 
mean-square pitch-angle deviations for co- and counter- 
propagating ensembles on these longer times would be 
misleading. Instead, as described above, we compare the 


values of the first minima of (or rather the slopes of 
the dotted lines through these minima). 

These values of A^^ for 21 different initial pitch-angles 
(/r(0) G {—1, —0.9,..., +1.0}) are then used to determine 
a best-fit parameter set s^, s“} for the 

quasilinear prediction for given by equation ( ^ . The 
remaining parameters (/cmin, 7 ) were fixed at fneir ac¬ 
tual values before the optimization was initialized. Al¬ 
though equation (1^ was derived under the assumption of 
a highly simplihea slab spectrum with a constant power 
index, hardly what our simulations of realistic MHD ex¬ 
hibit in Figure it will become clear that the agreement 
of A^^ and the quasilinear predictions of is rather 
striking. The results from a Levenberg-Marquardt fit us¬ 
ing our measurements of A^^(/r) are shown in Figure]^ 

The fit-values of match the values of a used for 
the cross-helicity injection well in all cases. A decrease of 
the fit quality if one Elsasser energy is very small (that 
is, for a = 0.8) is visible but not surprising. The rel¬ 
ative turbulence amplitudes SB^^^^/Bq are close to the 
actual value SB/Bq = 0.1, usually lying slightly below 
because the particle scattering is mainly due to only the 
slab component of the turbulence. The spectral indices 
obtained from the ht procedure are almost always larger 
than the predictions for the inertial-range spectrum in 
Kolmogorov’s or Kraichnan’s theories (5/3 or 3/2, re¬ 
spectively) since, like realistic turbulence, the spectra in 
our simulations steepen at the transition to the dissipa- 









Fig. 8.— Fully electrodynamic case: histograms of the pitch- 
angle cosine after ten gyroperiods for test-particles with t;(0) = 
and /i(0) = dil.O in a) balanced turbulence (cjc = 0.0), b) strongly 
imbalanced turbulence (ctc = 0.9) 

tion range (see Figure]^. 

The graphs of the quasilinear for these param¬ 
eter sets perfectly match the scattering rates de¬ 
termined from the test-particle trajectories, as shown in 
Figure Whereas the graphs in the balanced case are 
symmetric with respect to the sign of /r(0), particles in 
the imbalanced cases are scattered faster within the first 
gyroperiod if they travel in the direction of the weaker 
Alfven-wave population (if /r(0) > 0 for cr > 0 and vice 
versa). 

This asymmetry is shown more directly in Figure 
which contains pitch-angle histograms of particles with 
i'(O) = 3va after ten gyroperiods for test-particles ini¬ 
tially propagating exactly parallel or anti-parallel to the 
mean-field direction. In this case, both the heating 
that was observed for the negative-/i ensemble in Fig¬ 
ure and the initial magnetic moment are insignifi¬ 
cantly small, and the behavior predicted by quasilinear 
theory can be observed for much longer. While both 
ensembles, with ^(0) = -1-1 and with ^(0) = —1, are 
scattered equally fast in balanced MHD turbulence, pos¬ 
itive cross helicity leads to particles propagating along 
the mean-field direction being scattered much faster than 
in the opposite case. 

In Figure!^ we compare similar histograms of sample 
pitch-angle distributions after one gyroperiod with solu¬ 
tions that we directly obtained from the Fokker-Planck 
equation Q, using /(/i, t = 0) = i5(/i(0) — fi) as initial 
condition and assuming a homogeneous density along the 
z-coordinate. While we showed above that the mean- 
square displacement of the pitch-angle cosine after one 
gyroperiod agrees with the scaling of predicted 

by QLT, Figurej^demonstrates that the shape of the his¬ 
tograms agrees with a diffusive broadening of the pitch- 
angle distribution within the first gyroperiod. This con¬ 
firms once more that the quasilinear prediction ([^ gives 
an accurate description of the pitch-angle scattering on 
short timescales even in evolving imbalanced turbulence 
with electric fields. 

4.2. Magnetodynamic runs 

There are two possible causes for this asymmetry in im¬ 
balanced turbulence: the coherent interaction with prop¬ 
agating Alfven waves and the spatial structure of the 
turbulent fields. To identify the more important mecha¬ 
nism, we first ignore the electric-field component of the 
Lorentz force and propagate the test-particles using only 

Vi = avi X b(xi,t). (14) 


Fig. 9.— Fully electrodynamic case: histograms of the pitch- 
angle cosine after one gyroperiod for test-particles with a) t;(0) = 
VA and ^(0) = 0.0 in balanced turbulence with Uc = 0.0, b) i’(O) = 
3va and ^(0) = -1-0.5 in strongly imbalanced turbulence with ctc = 
0.9 
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Fig. 10.— Magnetodynamic case: evolution of Mmag and 
as in figure]^ 
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Fig. 11.— Magnetodynamic case: best-fit values for equation 
for for particles ensembles with t;(0) = va (left column) and 

u(0) = 3uy, (right), for four different values of a 

In the ‘magnetodynamic’ runs presented in this sub¬ 
section, the MHD fields were still dynamically evolved in 
parallel with the test-particles. Since the acceleration is 
perpendicular to the velocity at each time step, the ki¬ 
netic energy of each test-particle is perfectly conserved. 
As shown in Figure [TOl the evolution of both ensembles in 
the balanced case and of the positive-/! ensemble in the 
imbalanced case is similar to the fully electrodynamic 
setting presented above. 

However, the averaged magnetic moment of the en¬ 
semble with /!(0) = —0.6 visibly increases on longer 
timescales. Similarly, the distribution of M^ag contin¬ 
ues to broaden for at least ten gyroperiods, as indicated 
by the red shaded region in Figure [TOj:. Although the en¬ 
semble in Landau resonance with the dominant Alfven- 
wave population cannot be decelerated anymore because 
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H=COSV^/v H=COSV^/v 


Fig. 12.— Magnetodynamic case: pitch-angle scattering rate for cosmic-ray particles determined numerically after one gyroperiod 

and from a Levenberg-Marquardt fit to eqn. 0 See caption of figure]^ 


electric acceleration is not considered, the coherent inter¬ 
action with the magnetic component of the Alfven waves 
is now able to violate the approximate conservation of 
the magnetic moment. Energy is transfered from the Vz- 
component into the perpendicular velocity components. 
Without a constant adiabatic invariant, the pitch-angle 
cosine can decrease as before, even though the kinetic 
energy is fixed, and keeps increasing at a similar 

rate as in the fully electrodynamic runs. 

We use the first minima of A/i^ for p,(0) G 
{ —1, —0.9,..., -1-1.0} to determine a best-fit parameter 
set with the quasilinear in magnetodynamic turbu¬ 
lence for the same variables as before, but with given 
by equation Q. Most of the results are similar to those 
presented in me previous subsection, so we o nly show 
the best-fit cross-helicity values in Figure |ll| The 
results for the Alfven-speed ensemble (w(0) = va) are 
highly inaccurate, while the ensemble with ri(0) = 3va 
at least exhibits the correct trend: cr*-*^*^ grows weakly as 
a increases. 

Plotting and A^^ in Figure [l^ we see that, with¬ 
out electric fields, the scattering asymmetry is strongly 
reduced compared to the graphs in Figure 0 Electric 
fields obviously play a crucial role in establishing this 
asymmetry in imbalanced turbulence. 

4.3. Magnetostatic runs 

To determine the relative importance of wave propa¬ 
gation, we have also performed test-particle simulations 
with single static snapshots of the magnetic-field config¬ 
urations. In these magnetostatic runs, we calculate the 
Lorentz force on the particles using only the magnetic 
field, similar to the previous subsection (equation ([l4|)) 
but without evolving the field in time. 
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Fig. 13.— Magnetostatic case: Evolution of M^ag and as 
in figure]^ 


Since particles can no longer interact coherently with 
traveling waves, Landau resonance as observed before is 
impossible and pitch-angle scattering proceeds at simi¬ 
lar rates for positive and negative /t(0) even in stro ngly 
cross-helical turbulence. As shown in Figure the 
mean magnetic moment stays constant throughout the 
simulations and A/r^ (t) settles on an approximately con¬ 
stant value after a few gyroperiods. 

Evolving the trajectories for longer times, which is pos¬ 
sible since magnetostatic test-particle simulations are not 
very demanding compared to fully electrodynamic runs, 
yields a subdiffusive behavior of the pitch-angle distri¬ 
bution. Thus, we computed A^^ after one gyroperiod 
again to perform a best-fit operation, varying SB and 
s in equation ([^. As we found statistical noise to be 
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Fig. 14.— Magnetostatic case: pitch-angle scattering rate A^^(/i.) for cosmic-ray particles determined numerically after one gyroperiod 
and from a Levenberg-Marquardt fit to eqn. § See caption of figure]^ 


significantly larger in static runs than in the evolving- 
MHD runs presented above, we show a 10 %-error bar 
in the plots of A^^(/i). Repeating each test-particle run 
in ten different static snapshots of the magnetic field for 
otherwise identical parameters, we found that the values 
we obtained for varied by this amount, whereas the 
evolving-MHD runs only varied by about 2 %. 

Taking this statistical noise into account, the graphs 
of A^^(/i) are symmetrical with respect to the sign of 
/r (see Figure 141. Cross helicity is therefore irrelevant 
in magnetostatic turbulence, and the spatial structure 
of the magnetic fields can be ruled out as a contributor 
to the scattering asymmetry. As a sidenote, we remark 
that none of our results were affected by the sign of the 
test-particle charge, proving that magnetic-helicity con¬ 
tributions to the scattering asymmetry, which would be 
sensitive to the particle charge, are insignificant. 


4.4. Static-electromagnetic runs 

Adding the electric-field component back into the 
Lorentz force while still using only static field snapshots, 
we foun d that this asymmetry can be re-established (Fig¬ 
ure!^. While the absolute scattering rates A^^ are re- 
ducedTo about 70 % compared to the values for the fully 
electrodynamic runs in Figurej^ the shapes of the graphs 
are almost identical to their fully electrodynamic equiv¬ 
alents in each case, especially if one takes the statistical 
noise into account. 

Since the electromagnetic fields of the Alfven waves 
are fixed in time and space in these runs, it is the spa¬ 
tial structure of the electric fields that must be the pre¬ 
dominant soure of the scattering asymmetry. Hence, we 
conclude that particles in fully evolving MHD turbu¬ 
lence with non-zero cross helicity are scattered differently 


for positive and negative pitch-angles not only because 
of coherent wave-particle interactions, but additionally 
because the spatial structure of electric fields in imbal¬ 
anced turbulence favors such an asymmetry. Of course, 
this structure can be viewed as the result of the im¬ 
balance between Alfven-waves propagating parallel and 
anti-parallel to the magnetic mean-field direction, a prop¬ 
erty that is conserved when the fields are frozen. 

5. CONCLUSIONS 

We have investigated pitch-angle scattering of charged 
particles at velocities comparable to the Alfven speed in 
realistically evolving magnetohydrodynamic turbulence 
with non-zero cross helicity. The dependence of A^^ on 
the initial pitch-angle is qualitatively well described by 
quasilinear diffusion theory under the assumption of a 
simple slab spectrum. The agreement with the theoret¬ 
ical prediction is particularly impressive in the case of 
fully electrodynamic turbulence, showing that, for exam¬ 
ple, the interaction of Alfven-speed dust particles with 
shear-Alfven waves they encounter in the solar wind can 
be modeled accurately by equation ([5), at least on short 
timescales. An almost perfect fit of me quasilinear slab 
model to the numerical data is even possible in strongly 
cross-helical turbulence, as we have shown; therefore the 
scattering of charged particles in patches of strongly im¬ 
balanced turbulence, such as in the fast solar wind, is 
also described accurately by quasilinear theory. 

However, the degree of imbalance that we have to as¬ 
sume in order to obtain this perfect fit is lower than the 
imbalance which a direct analysis of the correlation of the 
turbulent MHD fields yields. Our results thus imply that 
quasilinear theory tends to overestimate the asymmetry 
of pitch-angle scattering in turbulence with extremely 
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high cross helicity. 

As we noted in Section pitch-angle scattering over 
small timescales, on whicn we have focused in this ar¬ 
ticle, is important for theories describing the injection 
of relatively low-energy ions into the Fer mi-acceleration 
mechanis m at quasiperpendicular shocks (Kirk & Heav- 
ens|1989 ). The structure of the shock-foot region is par¬ 
ticularly sensitive to the short-time behaviour reflected 
in the shape of A^^(/r). Hence, our results will bene¬ 
fit future theoretical work on quasiperpendicular shocks. 
Although in such applications the overall magnetic tur¬ 
bulence will be stronger than the intermediate levels that 
we have considered here, these results can easily be trans¬ 
ferred to pitch-angle scattering with regard to a magnetic 
mean-field averaged only over intermediate lengthscales. 
On such smaller scales, the aniso tropy with respect to the 
mean-field is more pronounced ( Howes et al.||20lT ), and 
our results are valid. Eventually, the local pitch-angle 
diffusion thus derived could be combined with models 
for fieldline diffusion and crossfield diffusion to obtain a 
self-consistent description of the shock-f oot region, con¬ 
ceptu ally similar to compound diffusion (Kota & Jokipii 
200 ^. 

IFshould be noted that our simulations have focused 
on incompressible magnetohydrodynamics and did not 
include transit-time damping by the fast magnetosonic 
wave. The compressional fast mode has been shown to be 
more efficient at accelerating and scattering charged par- 


ticles (Schlickeiser & Miller 1998 Yan fc Lazarian|[2602 1 
than the shear-Altven modes included in this investiga- 
tion. However, the assumption of incompressible MHD 
and hence of a negligible contribution of the fast mode 
is still valid in many astrophysical systems, such as the 
fast solar wind. 


In our magnetodynamic simulations, which differ from 
the previously described fully electrodynamic case in that 
we have neglected the acceleration force due to the elec¬ 
tric fields, the pitch-angle asymmetry of the scattering 
during the hrst gyroperiod is significantly smaller. In 
static-electromagnetic runs, in which the electric force on 
the test-particles was included but the turbulent fields 
were frozen in time, this asymmetry is almost as pro¬ 
nounced as in the electrodynamic runs. This contrast 
implies that the pitch-angle dependence of the scatter¬ 
ing coefficients in imbalanced turbulence is mainly due 
to the spatial structure of the electric fields, and only 
secondarily due to coherent wave-particle interactions, 
with possibly important consequences for the isotropiza- 
tion rate of cosmic-ray particles in shock fronts. It is 
clear from these observations that the predictive power 
of quasilinear theory for the diffusion of cosmic-ray par¬ 
ticles in imbalanced turbulence profits immensely from 
using a more realistic spectral description of cross-helical 
fields. 

On timescales exceeding two or three gyroperiods, we 
observed that the mean-square deviation of the pitch- 
angle cosine settled on an almost constant value for most 
test-particle ensembles, in contradiction to the diffusive 
broadening one would expect. This spreading continues 
for significantly longer if the test-particles are close to 
Landau resonance with the dominating Alfven-wave pop¬ 
ulation in imbalanced turbulence. While the presence of 
electric fields allows for an energy change that accom¬ 
panies the longer-lasting spread of A/i^, we found that 
magnetodynamic runs conserve energy perfectly, but vi¬ 
olate the adiabatic invariance of the magnetic moment 
for such test-particles. Hence, without accou nting for 
adiabatic focusing as described, for example, in Schlick- 
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eiser & Shalchi 

(2008) and 

Schlickeiser & Jenko ( 

2010 


the quasiiinear description ot pitch-angle scattering in 
imbalanced turbulence is valid, but only for a few gy- 
roperiods. 
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APPENDIX 

THE ONE-PERIOD RESONANCE FUNCTION TliTc) 

In standard applications of quasiiinear theory, the resonance condition ujr = ± is obtained from the infinite¬ 

time limit of the resonance function TZ{t). For example, the pitch-angle diffusion coefficient in the magnetostatic- 
turbulence model is 


where the resonance function is defined as 


= lim / dfc^ 




R2 

-^0 


7Z{t) = 3? / drexp(izz7T — Pr), 

Jo 


(A2) 


with w = (jjR — kzVz^^- Here F denotes the damping rate of the oscillations. For general t, performing the integration 
and taking the real part yields 

^ F [1 - e-^ * cos(CTt)] +w sm{wt) 


F^ -I- 

the infinite-time limit of which (Ft ^ 1) we denote as TZac,: 


F 


7^oo = lim 7^(^) = 2 , 2 • 

t—^oo i ^ -p 

Thus, the resonance function can be written as 

7?,(t) = 77,00 + 77.corr(0: 

where we have defined a correction term as 

„ , , _rt 137 sin(ci7t) — F cos{'cut) 

'VcorrUj = e TOO , 9 ■ 

In the weak-damping limit (F —> 0), TZao reduces to a Dirac function, 

lim 77-00 = TT 6(w), 
r-j-o 

such that performing the integral in equation (Al) yields the standard result stated in equation (H). For finite times, 
however, 77,corr(Z) gives an additional contribution even in the limit F —>■ 0: 


(A4) 

(A5) 

(A6) 

(A7) 


.TO / N N sm(ci7t) 

lim TZ{t) = TT S{w) -I-. 

r-fO va 


(A8) 


Since the latter term is positive for t = 2b, the quasiiinear result 27^^ represents a lower limit for the pitch-angle 
scattering rate defined in equation (10). 
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